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Division of Applied Mathematics 
Brown University 
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ABSTRACT 

In this paper we discuss the wave-resolution properties of the Fourier approximations of 
a wave function with discontinuities. It is well known that a minimum of two points per 
wave is needed to resolve a periodic wave function using Fourier expansions. For Ohebyshev 
approximations of a. wave function, a minimum of 7T points per wave is needed [3]. Here we 
obtain an estimate for the minimum number of points per wave to resolve a discontinuous 
wave based on its Fourier coefficients. 

In our recent work on overcoming the Gibbs phenomenon, we have shown that the Fourier 
coefficients of a discontinuous function contain enough information to reconstruct with expo- 
nential accuracy the coefficients of a rapidly converging Gegenbauer expansion. We therefore 
study the resolution properties of a Gegenbauer expansion where both the number of terms 
and the order increase. 


1 Research supported by AFOSR grant 90-0093, ARO grant DAAL03-91-G-0 123, NASA grant NAG 1-1 145, 
DARPA grant N00014-91-J-4016, and by NASA contract NAS1-18605 while the authors were in residence at 
ICASE, NASA Langley Research Center, Hampton, VA 23665. Computation supported by the Pittsburgh 
Supercomputer Center. 


1 




1 Introduction 


In [5] (see also [3]), the issue of the number of points required to resolve a wave has been 
considered as a measure of the accuracy of a given scheme. In particular it has been shown 
that only two points per wave are required if the wave function is approximated by a Fourier 
series. This result is almost trivial since the expansion basis contains the wave function. It 
is also shown in [3] that if the wave function f(x) = e l ™ x is approximated by a truncated 
Chebyshev expansion then a minimum number of tt points (or, alternatively, number of terms 
retained in the expansion) per wave is required. The proof is based on the observation that 
if the number of terms in the expansion exceeds this minimum of 7 r per wave then the error 
decays exponentially (see [3, page 35]). 

The situation changes when the wave function is discontinuous. As a generic example, 
consider the above mentioned function f(x) = e t7TuJX defined on [— 1 , 1 ] where u> is not an in- 
teger. In this case, the truncated Fourier expansion does not converge at all in the maximum 
norm. This is known as the Gibbs phenomenon and there is no meaning to the question of 
number of points required to resolve such a wave. 

Recently, however, the Gibbs phenomenon has been resolved. In [4] it has been shown 
that the first N Fourier coefficients of an analytic, but nonperiodic function contain enough 
information on the solution to construct an exponentially convergent Gegenbauer series. 

Consequently we are concerned with resolution properties of the expansion in Gegenbauer 
polynomials C^(x) to approximate a non-periodic wave function f(x) — e l7VUJX where u; is not 
an integer. We first consider the question of minimum number of points per wave needed 
to resolve this function. Two separate situations are considered. The first involves a fixed 
A in the Gegenbauer expansion. Note that the Chebyshev expansion is a special case with 
A = 0. For this situation our results show that the minimum number of points per wave is 
7T, the same as in the special Chebyshev case. In the second situation we assume that 7 = ~ 
is a constant in the Gegenbauer expansion. We show that the minimum number of points 
per wave increases with 7 at the boundary x = ±1 (i.e. at the discontinuity of the wave) 
and decreases with 7 at the center x = 0. This second situation corresponds to the case in 
which the first N Fourier coefficients of the function are known, since in this case we can 
reconstruct the Gegenbauer expansion for 7 = — constant. 

In Section 2 we quote some results about Gegenbauer polynomials. In Section 3 we study 
the resolution properties of the Gegenbauer expansion of a wave function, first for A fixed 
and then for 7 = ~ constant. Section 4 is devoted to the discussion of how many Gegenbauer 
coefficients can be obtained accurately from the first N Fourier coefficients. Here we have 
the minimum number of points per (discontinuous) wave to get exponential decay of the 
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error. Section 5 brings numerical confirmation for the above results. 

We will use A to denote a generic constant independent of all the growing parameters 
throughout this paper. The actual value of A may be different at different locations. 


2 Preliminaries 


In this section we will collect some definitions, equalities and inequalities which will be used 
in later sections. 

We start by defining the Gegenbauer polynomials C^(x). 

Definition 2.1 For A > 0, the Gegenbauer polynomial C*(x) is the polynomial of degree n 
that satisfies 

J (1 — x 2 ) x ~^Cl(x)C^{x)dx = 0 k / n (2.1) 

with 


<£( 1) = 


T(n + 2A) 
n!T(2A) 


( 2 . 2 ) 


The Gegenbauer polynomial C^(x) achieves its maximum at the boundary x = 1 [2, page 


206] 


\C x n {x)\ < C x n { 1) 


and its value at the center x = 0 is given by [1, page 777] 


(2.3) 


c£,(0) = (-i) 


„ T(A + n) 
n!T(A) ’ 


c 2 A ntl (0) = 0. 


(2.4) 


The Gegenbauer polynomials thus defined are not orthonormal. The norm of C*(x) is 
given by [2, page 174] 


K = j\ i - x^-ic>(x)ci(x)ix = 

We will also need the following identity, which can be found in [2, page 213] 


(2.5) 


_L i _ x 2 ) x ~h'™*C?{x)dx = T(A) (“) A + A) J/ + a(ttu;). (2.6) 

Throughout this paper we repeatedly use the Stirling’s formula: for any x > 1, 

T(x + l) < {2n)^x x+X e- x e^ 
r(z + l) > (2ir)iat' + ie _ * . 
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(2.7) 

( 2 . 8 ) 




3 Wave Resolution Properties of Gegenbauer Expan- 
sions 

Consider the non-periodic wave function 

/(x) = e* w (3.1) 

where u> > 0 is not an integer. The Gegenbauer coefficients of this function are defined by 

!\ 0 = Yx j\ 1 - x 2 ) x -te'™ x Cf(x)dx (3.2) 

with hf given by (2.5). 

Our objective is to find a lower bound on the ratio 

m 

r = — (3.3) 

where u is the number of waves in (3.1) and m is the number of terms in the Gegenbauer 
expansion 
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(3.4) 


/i(*) = E/‘(OCi A (*) 

/=0 

such that the approximation error using the expansion (3.4) is exponentially small when 
m -* oo. The ratio r defined in (3.3) is usually called the number of points per wave. It is 
literally the number of points per wave for the collocation case and can be called number of 
modes per wave for the Galerkin expansion (3.4)-(3.2). We can then define the regularization 
error RE(\,m,r,x ) to be 

RE{ A, m, r, x) = | /(*) - &{x ) | (3-5) 

and ask the question of finding a lower bound r 0 = r 0 (A,x) such that the regularization error 

(3.5) is exponentially small for r > r 0 when m — * oo. 

According to the identity (2.6) and the definitions (3.2)-(3.4), the regularization error 

(3.5) can be explicitly expressed as 

RE( A,m,r,*) = r(A) (— V £ i'(l + X)J M (^) Cftz) . (3.6) 

l=m + 1 V r ' 

first we have the following theorem for the case with fixed A. 

Theorem 3.1 If A is fixed, then r 0 = 7T is a lower bound for the number of points per wave 
to obtain exponentially small regularization error for all -1 < x < 1 . 

Proof: Assume r > tt. Denote z, = and z = f. For Z > m we have z, < 2 < 1 and 

hence q(z,) < q(z) < 1 according to Lemma 2.1. Take 6 = ±±2^ < 1 and m big enough so 
that ( '|/ + + A A ){i+ 4 i r if z( z ) < 5 for l > m. If we define 

B(l) = (( + A) (,(,))'“ (3.7) 

we have, for l > m 

g ( /+1 ) = ( /+1 + A )( /+ ^l g (z) < s < 1 . (3.8) 

B(l) (Z + A)(Z + 1) 

We can then start from the explicit formula (3.6) and absorb all the A dependent terms 
into the generic constant A to obtain 

RE(\,m,r,x) < Am ^ (l + A) Ji+\ ( — — J C t (1) 

/=m+l V ' 

< Am'* £ (! + A)«*))' +aJ T±^> 
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< 

< 

< 


oo 

Am~ x J2 B ( l ) 

[=zm+\ 


Am 


-a B{m) 


Am' (r/(z)) 


1 — s 

m+A 


where we have used (2.3) in the first inequality (2.9) and (2.2) in the second inequality, the 
definition of B(l) in (3.7), the fact that q(z{) < q(z) in the third inequality, and (3.8) in the 
fourth inequality. This finishes the proof because q(z) <1. □ 


The result of Theorem 3.1 is not surprising; for A = 0 this is simply the known result for 
the Chebyshev case [3]. 

For fixed A, the estimates are essentially the same for the boundary point x = 1 or for 
the center point x = 0, since there is only an algebraic difference between (7^(1) and C^O). 
However, this algebraic difference becomes bigger when A increases. 

For our purpose we are more interested in the case A ~ rn. This is because we are 
interested in the situation where the finite Fourier series of f(x) is given. In [4] we proved 
that it is possible to recover uniform exponential accuracy from the finite Fourier series of a 
non-periodic analytic function through the use of Gegenbauer expansions with A ~ m. 

Let us now assume that 7 — ~ is a constant. We have then the following theorem. 
Theorem 3.2 If 7 = — is a constant, then 

' m “ 

(1) The regularization error RE^m, rn, r, 1), at the boundary x = 1, is exponentially 
small if r > max(rq,r 2 ) where 77 is the unique solution of 


P 


7T 


( 2 « 7 ( 1 + 7 ))^ |?i(t) 


in the region 77 > 


r i(l +7)/ (1 + 27) 1+2 ^ 

if (3.9) has a solution (i.e., if p(c( 7)) > ^1(7)), or r i 


(3.9) 


n/T+27 


otherwise. Also, is the unique solution of 


7 r 


72(1 + 7)/ I+27’ 


(3.10) 


(2) The regularization error RE^m, m, r, 0), at the center x = 0, is exponentially small 
if r > max(r 3 ,r 4 ) where r 3 is the unique solution of 


P 


( * 

v 3(1 +7) 


(e(l + 7))^ 

(1 + 27) 2 +7 


= QM 


(3.11) 


within the region r 3 > if (3.11) has a solution (i.e., if ^(0(7)) > Q 3 ( 7)), or r 3 = 

otherwise. Here r 4 is the unique solution of 
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7T 


1 


HmIT T)) = 7tt^- (312) 

The functions p(z ) and q(z) in (3.9)-(3.12) are defined by (2.10) in Lemma 2.1 and the 
constant 0(7) is defined by (2.11). 

Proof: Denote z\ = an< ^ z ~ r(m+ A) = r (i + 7 ) • For l > m we have zt < z < c( 7) when 

r > r i > hence 7(2/) < 7(2) < according to Lemma 2.1 and the condition (3.10). 

Take s = i±lL±|2Mfl ) which yields (1 + 27)7(2) < s < 1. Hence we can take m big enough 
so that for / > m, 


(/ + 1 + 7m)(i 4 - 27 m) 
(l -f 7 m)(l + 1) 


q(z) = (l + 


1 


/ \ ( l + 27m' 


< 


l + 7m J \l + 1/ \ / + 27/ 

1 + 77^) (1 + S * < 1 • 


(1 + 27)7(2) 


If we now define 


B(l) = (l + im ) (<,(z))' +TO , 

we have for / > m 


B(/+l) (/+ 1 +7m)(/ + 27m) 

~mr = </+ 7 m)(i+i) ,(2)£j<1 

We can then start from the explicit formula (3.6) to obtain 


(3.13) 


(3.14) 


RE(ym, m, r, 1) < 
< 
< 
< 
< 


r(7 m ) ( 


2 r \ 
\ 7 rmJ 


) S + 7 m ) 


/rrm+l 


■w (™) 


(7^(1) 


r( 7 m) f2ry m ^ ni w , w+7m r(/ + 2 7 m) 

mim) Vim) E +i (' + 7m)(,(z,)) — jj— 

r( 7 m) / 2r y- - 

T(27m) \7rm/ 
r( 7 m) / 2r y m g(m) 

T(27m) \7rm/ 1 

< t (1 +2 7 ) ,+i ’ 


s 

v m 


,(2c7(l +7))' 


p(z) 


where we have used (2.9) and (2.2) in the second inequality, the definition (3.13) and the 
fact 7(2/) < 7(2) in the third inequality, (3.14) in the fourth inequality and Stirling’s formula 
(2.7)-(2.8) in the last inequality. 
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According to condition (3.9) and Lemma 2.1, the number ^(^+ 7 ))^ ( r ( i+-y) ) strictly 
less than one for r > 77. This finishes the proof for x = 1. 

The proof for the center point x — 0 is similar and is thus omitted. □ 




Figure 1 : Minimum number of points per wave as a function of 7 . Left: minimum — ~ 
to obtain exponentially small regularization error; Right: minimum R = ~~ to obtain expo- 
nentially small resolution error from Fourier partial sum. 

In Figure 1 (left) we show the curves of number of points per wave, r = versus 7 , 
both for resolving the boundary x = 1 (upper line) and for resolving the center x = 0 (lower 
line). These were obtained from the results (3.9)-(3.12) in Theorem 3 . 2 . We can see that 
the top curve for x — 1 is an increasing function of 7 , while the bottom curve for x = 0 is a 
decreasing function of 7 for 7 > 0.37. At 7 = 1 , we need approximately 7.03 points per wave 
to resolve the boundary and 3.51 points per wave to resolve the center. At 7 = 5 , this two 
numbers change to 8.11 and 2.33, respectively. Figure 1 (right) is related to the truncation 
error to be discussed in next section. 


4 The Truncation Error from Finite Fourier Series 

In this section we consider the situation that f(x) = | s not known and only the first 

2N + 1 Fourier coefficients f(k) are given. We would like to recover, within exponential 
accuracy at all points in — 1 < x < 1 , the finite Gegenbauer expansion f^(x ) of f(x) as 
defined by (3.4). This, together with the results in the previous section, will establish the 
number of points per wave for resolving discontinuous waves based on its Fourier partial sum 
through Gegenbauer polynomials. The result in this section parallels the general theorems 
in [4], 


7 





The Fourier partial sum of the discontinuous wave function f(x ) = e tnwx is given by 

M*) = E /(*)«"” (4.1) 

k=-N 

where the Fourier coefficients f{k ) are defined by 

(4.2) 

Assume that the first 2 N + 1 Fourier coefficients f(k), |fc| < N and hence also the 
Fourier partial sum (4.1) are known but the function f(x) is not. We thus do not know the 
exact Gegenbauer coefficients f x (l ) of f(x) defined by (3.2), but only the approximate ones 
obtained from the Fourier partial sum /n(x), which we denote by g x (l) 


9 X {1) = - x2 ) A *fN(x)C?(x)dx. (4.3) 

Notice that g x (l ) depends on N. 

We now define the truncation error T N,x) to be 


T E( A, m, N, x) = 


E(f x 0) - ^(0)0^) . 


1=0 


(4.4) 


Notice that u> or r = ~ is not an explicit parameter in our definition of the truncation 
error. The estimates we obtain later will be uniformly valid for all u>. 

In the next theorem we shall bound the truncation error in terms of N, the number of 
given Fourier coefficients, m the number of terms in the Gegenbauer expansion, and A. 
Theorem 4.1 The truncation error at the boundary x = 1 satisfies the estimate 


T E(X, m, N, 1 ) < A 


(m + A)T(m + 2 A)T(A) ( 2 


(m - l)!r( 2 A) 

At the center x = 0, the truncation error satisfies 


7T 


N 


A-l 


(4.5) 


TE(X,m,N,0) < A 


(ro + A)r(f + A) / 2 


r (f) 


(-) 

\ttNJ 


A-l 


(4.6) 


If A = 7 m and m = /3N where 7 and /? are positive constants, then the truncation error 
satisfies 


T E{~ifiN,(iN, N, 1) < AN 2 


7^(1 + 2 7 ) 1+2V 
, (27re7)'> 


0N 


at the boundary x = 1 and 


(4.7) 
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TE^N, /37V, TV, 0) < AN 2 (4.8) 

at the center x = 0. 

Proof: Si nee 

f(x) - f N (x) = £ f(n)e' n ” X (4.9) 

|n|>/V 

and \f{n)\ < 1 according to the definition (4.2) with f(x) = e t7ru;:r , we can estimate the 
truncation error as follows 


TE(\, m, N, x) < (m + l) o max (/ A (/) - ^ A (/))C' / A (x)| 


(m + 1) max 

' 0 <l<m 


= (m+ 1) max 


o</< 


£ hn)C?(x)[-J (1 — x 2 ) x ~ie iniry C^(y)dy 

n | >N \ n l - 7-1 

E /(")C, A (x) r(A) ( 

,, I ^ AT ' 


|n|>7V 

2 ' A_1 


— - ) i‘(l + \)Ji+\(Tm) 

s 7rn/ 


< /lmr(A)(— ) max(/ + A)|(7 A (x) . 


(4.10) 


In the second step we have used (4.9), the definition of f x (l) in (3.2) and that of <) A (/) in 
(4.3), in the third step we have used the equality (2.6), and in the last step we have used 
the facts |/(n)| < 1 and |J„(x)| < 1 for all x and u > 0 [1, page 362]. 

For the boundary point x = 1 we can then proceed as follows 


/ 9 \ A_1 

TE(\,m, N, 1) < /bnr(A) ( -=-) 

KnN ) 




= Am 


r(A) / 2 \ 


A— 1 


max 

0</<m 


(/ + A)r(/ + 2 A) 


< 


r(2A) Vt tNJ 
(m + A)r(m + 2A)r(A) / 2 \ 
(m — l)!r(2A) \Vn) 


l\ 

A— 1 


(4.11) 


where in the second step we have used the formula (2.2) for C A (1), and in the last step 
we have used the fact that ( /+A ) I 4 ;+2A ) an increasing function of l. The result in (4.5) is 
thus proven. Some simple algebra and the use of Stirling’s formula (2.7)-(2.8) easily produce 
(4.7). The estimates (4.8) and (4.10) at the center point x = 0 can be obtained in a similarly 
fashion and the detail is thus omitted. □ 


Since our estimates of the truncation error, (4.7) and (4.8), do not depend on the wave 
number u or r = the minimum number of points per wave for exponential convergence, 
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defined by R = is obtained by using the r = ^ derived in Theorem 3.2 and the largest 

/? in (4.7) or (4.8) such that the factors on the right hand sides which are being raised to 
power N are still less than one. This gives the following theorem. 

Theorem 4.2 Assume that the Fourier coefficients 

f(k) = j /_) f(z)e-“"*dx 

of the function f(x) = e wirx are known for — N < k < N. Let g x (l), 0 < / < m be the 
Gegenbauer expansion coefficients of fiv{x) = Ylk=-N f{k)e' k * x given by (4.3). Define the 
number of points per wave by 


and the resolution error by 



c J 


(4.12) 


E( A, m, N, R , x ) = 


m 


(=0 


Then, if A = 7 m and m = f3N, we have the following results. 
(1) At the boundary x = 1, if one uses 


(4.13) 


0< 


27re7 

(1 + 2 7 )^ +2 


= B,( 7 ) 


(4.14) 


and r = ^ > max(ri,r 2 ) as is derived in Theorem 3.2, then the resolution error is exponen- 
tially small and bounded as 


E(y/3N, 0N, N, R, 1 ) < A (N 2 b% + 1 $) 


(4.15) 


where 


b? — 


'/T(l +2 t) 1+27> 


<i, ^ 


■( 


(1 + 2 7 ) 1+2 ' y / 7T 


\ (27re 7 )' Y 

The number of points per wave 7? = ^ is estimated by 

2max(ri,r 2 ) 


(2e 7 (l + 7 )) 7 \r(l + 7 ) ; 


/? > 


Bi(7) 


where 77, r 2 are defined in Theorem 3.2 and 7?i( 7 ) is defined by (4.14). 
(2) At the center x = 0, if one uses 


P< 




(1 + 2 7 ) 2 ^ 


ET7 = 


< 1 . 


(4.16) 


(4.17) 
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and r — ^ > max(r 3 ,r 4 ) as is derived in Theorem 3.2, then the resolution error is exponen- 
tially small 


E(yftN, ftN, N, R, 0) < A (N 2 c% + <%) 


where 


ct = 


>(1 +2 7 )5+ 7 ' 


(1 + 27)2 +7 


The number of points per wave R = — is estimated by 


R> 


2max(r 3 ,r 4 ) 

7 ) 


where r 3 , r 4 are defined in Theorem 3.2 and # 2 ( 7 ) is defined by (4.17). 
Proof: This is simply a combination of Theorems 3.2 and 4.1. 


(4.18) 


(jre)’ I CR 1 (e(l + 7 ))’ P (r(l +7),^ ' 


(4.19) 


□ 


In Figure 1 (right) we show the curves of number of points per wave, R = versus 7, 
both for resolving the boundary x = 1 (upper line) and for resolving the center x = 0 (lower 
line), obtained from the results (4.16) and (4.19) in Theorem 4.2. We can see that the top 
curve for x = 1 achieves its minimum at around 7 = 0.9. and the bottom curve for x = 0 
is an increasing function of 7. We emphasize that these results are obtained with different 
ft for x = 1 and x = 0 (from (4.14) and (4.17) respectively). In practice a single ft should 
be used since one would like to resolve both the boundary and the center simultaneously. If 
the ft for x = 1, given by (4.14), is also used for x = 0, the minimum number of points per 
wave, 77, to resolve the center 1 = 0 , would be described by the middle curve in Figure 1 
(right). For the single ft chosen according to (4.14), at 7 = 1, a minimum of 22.2 points per 
wave is needed to resolve the boundary x — 1, and a minimum of 11.1 points per wave is 
needed to resolve the center x = 0. 


5 Numerical Results and Conclusions 

In this section we perform numerical calculations to demonstrate the theory developed in 
previous sections. We use the discontinuous wave function 

f(x) = cos(7ru;(£ -f 1 )) ( 5 . 1 ) 

for various wave numbers u and report both the regularization error defined by (3.5) and 
the resolution error defined by (4.13). 
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We implement the method in the following way. The exact Gegenbauer coefficients of 
f(x), which are needed for the regularization error, are computed using (2.6). The approx- 
imate Gegenbauer coefficients g x (l ), defined in (4.3), which are needed for the resolution 
error, are computed using the following formula 

g x (l) — <^o;/(0) + + A) £ J/ + a(ttA:) rd f(k) (5-2) 

0<\k\<N V7rfc/ 

where f(k ) are the Fourier coefficients of f(x) defined by (4.2). This formula can be easily 
derived from the definition of g x (l) in (4.3) and the integration formula (2.6). We compute 
the Bessel function J„(x) using an IMSL routine. The approximation to f(x) is obtained by 
directly summing (3.4) for the regularization error, or by directly summing 


<£(*) = £* A ( 0<?(*) (5-3) 

/=0 

for the resolution error. The Gegenbauer polynomials C^(x) are computed by the formula 


C/ A ( cos 


i 


*) = £ 

k = 0 


r (k + A) r{i - k + a) 
fcir(A) (/ — A.- ) ! r ( A ) 


cos(/ — 2 k)d 


(5.4) 


which can be found in [1, page 175]. 

We remark that the implementation techniques described above are subject to roundoff 
effects for large A and m. A better way of implementing the method through Chebyshev 
polynomials is currently under investigation. 

In Figure 2 we show the errors, in a logarithm scale, at the discontinuity x = 1 for u = 1.4, 
2 .4, 3.4 and 4.4. We choose 7 = ^ = 1 in this illustration. On the left, the regularization 
error RE(m,m, ^,1) is shown as a function of m, the number of terms retained in the 
Gengenbauer expansion. On the right, the resolution error £(y, y, Af, 1) is shown as a 
function of 2N, to total the number of terms in the Fourier expansion. Here we take ft = 0.5 
which satisfies (4.14). We can see that the errors are order 0(1) until r = ^ or fi = y 
reaches the critical values obtained in Theorem 3.2 and in Theorem 4.2, after which the 
errors drop exponentially. Figure 3 shows the same result but for the center point x = 0. 
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Figure 2: Logarithmic error at the discontinuity x = 1 for (5.1) with u = 1.4, 2.4, 3.4 and 
4.4. Here 7 = £ = 1. Left: the regularization error RE(m,m,^, 1) as a function of m; 
Right: the resolution error E(j, y, N, 1) as a function of 2 N, ft = 0.5. 



Figure 3: Same as Figure 2 but for the center point x = 0. 
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